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We recall results concerning one-dimensional classical and quantum systems 
with ladder operators. We obtain the most general one-dimensional classical 
o j systems respectively with a third and a fourth order ladder operators satis- 

fying polynomial Heisenberg algebras. These systems are written in terms of 
the solutions of quartic and quintic equations. They are the classical equiv- 
alent of quantum systems involving the fourth and the fifth Painleve tran- 
scendents We use these results to present two new families of superintegrable 
systems and examples of trajectories that are deformation of Lissajous's fig- 
ures. 



1 Introduction 

The most well known system with ladder operators is the ID harmonic os- 
cillator with first order raising and lowering operators. Such operators first 
order creation and annihilation operators were also used in context of the 
N-dimensional isotropic and the anisotropic harmonic oscillator. The inte- 
grals of motion of these mult i- dimensional systems can be written in terms of 
these operators [1,2,3]. The second order ladder operators of two-dimensional 
superintegrable Hamiltonians were obtained in the 60ties by Winternitz and 
al [4] and a classification of two-dimensional systems with second order lad- 
der operators was performed later by Miller and Boyer [5]. More recently, 
systems with third order ladder operators were studied [6-19]. It was shown 
that a ID quantum system with third order ladder operators can be con- 



1 



structed from supersymmetric quantum mechanics (SUSYQM) with a first 
and second order supercharges [11,16,19]. The potential, the supercharges, 
the ladder operators and the zero modes were written in terms of the fourth 
Painleve transcendent (P4). The case of quantum ID systems with fourth 
order ladder operators and the relation with the fifth Painleve transcendent 
was first explored in Ref.20 and Ref.21 and also later with extensions to 
2D superintegrable systems and the relations with integrals of motion [22]. 
These systems were written explicitly in terms of the fifth Painleve transcen- 
dent (-P5). These systems with higher order ladder operators can allow many 
zero modes for the annihilation and creation operators. It was discussed 
how superintegrable systems generated from such ID systems can have very 
interesting spectrum with more complicated degeneracies [16,22] due to the 
presence of singlet, doublet and triplet states. 

In classical mechanics, a particular case of a system with third order integral 
of motion [13,23] was related to third order ladder operators [24]. Let us 
mention that classical systems with third order ladder operators were also 
explored in an earlier paper of Eleonsky and Korolev [18]. The purpose of 
this paper is to present a classification of classical systems with third and 
fourth order ladder operators satisfying polynomial Heisenberg algebras and 
show how these results can be used to generate superintegrable systems. 

The polynomial (i.e. polynomial in the momenta) ladder operators de- 
scribed in Ref.11,16,18 and 19 satisfy polynomial Heisenberg algebras. How- 
ever, for many quantum systems we need to introduce nonpolynomial ladder 
operators and generalized Heisenberg algebras. For example, ID systems 
with quadratic energy spectrum and their classical equivalent allow such 
nonpolynomial ladder operators [18,25-28]. The infinite well and the Morse 
potentials are examples of such systems. The existence of raising and lower- 
ing operators is also not restricted to systems with separation of variables in 
Cartesian coordinates [3] and Hamiltonians with a scalar potential. Ladder 
operators in context of systems with a monopole interaction were discussed 
recently in the case of generalized MICZ-Kepler [29] and generalized Kaluza- 
Klein monopole [30]. 
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These ladder operators are also related to other approaches in quantum 
and classical mechanics. The relations between integrals, supersymmetric 
quantum mechanics and ladder operators were discussed by various authors 
[1-5,14,15,22,31-36]. These connections provide new methods to study know 
superintegrable systems, or even obtain new one and their corresponding 
polynomial algebra [14,15,22,36]. 

Let us present the organization of the paper. In Section 2, we recall results 
on ID systems with ladder operators in quantum and classical mechanics. In 
Section 3, we obtain classical systems with third order ladder operators that 
satisfy a polynomial Heisenberg algebra. In Section 4, we discuss the case 
of fourth order ladder operators. In Section 5, we recall a method discussed 
in earlier paper to generate superintegrable systems in context of classical 
mechanics from analog of lowering and raising operators [14]. We obtain 
from results of Section 3 and Section 4 two new families of superintegrable 
systems and moreover we present examples of trajectories. 



2 Quantum and classical systems with ladder 
operators 

Let us consider the following ID Hamiltonian of the form : 



H=^- + V(x), P=-ihd x , (2.1) 
with a polynomial ladder operators of order n (i.e. polynomial in momenta), 

.±iP n + f n _ 1 (x)P n - l + ...±if 2 (x)P 2 + f l (x)P±if even, 
^ P n ± 2/ n _ 1 (x)F™- 1 + ... ± if 2 (x)P 2 + h{x)P ± if odd, 



that satisfy a polynomial Heisenberg algebra in quantum mechanics or the 
classical equivalent Poisson algebra in the context of classical mechanics. 
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In quantum mechanics, a polynomial Heisenberg algebra satisfy the follow- 
ing algebraic relations : 

[H, A^} = uA\ [H, A] = -uA, (2.3) 

AA^ = Q(H + u), A^A = Q(H), [A,A*]=P(H), (2.4) 

with P{H) = Q{H + uS) — Q(H). The polynomial Heisenberg algebra gener- 
ated by {H, A, A^} provides informations on the spectrum of H. We have for 
Q{H) a n th-order polynomial in H . The annihilation operator can allow at 
most n zero modes (i.e. a state such Aif) = 0). In each axis we can have at 
most n infinite ladders by acting iteratively with the creation operator. The 
creation operator can also allow zero modes (i.e. a state such A'tp = 0) and 
we can have finite ladder operators. Such operators and their corresponding 
polynomial Heisenberg algebras are important in context of coherent states 
[37]. 

In classical mechanics we can also define analog of raising and lowering 
operators [2] (A~ and A + ) and the corresponding polynomial Heisenberg 
algebra. We replace commutation relations by Poisson brackets. These op- 
erators satisfy the algebraic relations : 

{H,A + } p = uA + , {H,A~} p = -uA- , (2.5) 

{A~, A + } p = P(H), A~A + = A + A~ = Q(H), (2.6) 

where P(H) and Q(H) are polynomials. They are related by P(H) = 
—iuj^^p-. Generalizations of the algebraic structure in classical and quan- 
tum were obtained by considering u = u(H). Ladder operators can also 
be useful in context of ID systems and they can be used to generate time- 
dependent integrals of motion : 

Q± = A ± e^\ (2.7) 

^ = WQ±}p+ ^ = . (2.8) 

The constant value of this integrals given by Eq.(2.7) can be written in the 
following form : 

q ± = c(E)e ±ie ° . (2.9) 
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We have A ± as given by Eq.(2.2) and thus 



A ± = c(E)e ±{9o+u)t \ (2.10) 

generate equations that can be used to obtain the trajectories on the phase 
space (x(t), P(t)) algebraically [27]. Ladder operators are thus useful in con- 
text of quantum mechanics but also classical mechanics. 



2.1 Classification of systems in classical and quantum 
mechanics with ladder operators 

The most general system with first and second order ladder operators in 
classical and quantum mechanics coincide [5] : 

First order : 




(2.11) 



Second order : 

V 2 = + |, (2.12) 

These potentials also appear in context of superintegrable systems [4]. Re- 
cently systems with third and fourth order ladder operators were obtained 
and studied [11,16,21,22]. 

Third order : 

V 3 (x) = ^x 2 + —P 4 + —P* + uVh^xP 4 + —(-a + e) , (2.13) 

with (P 4 = P 4 (yffx, a, 0) and e = ±1). 
Fourth order : 

KJ 8 V (P 5 -1) 2 P 5 ' x 2K 8 P 5 ' V 2(P 5 -1) 

(2.14) 
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with (P 5 = P 5 (j-x 2 , a, b, c, — |) and e = ±1). The second order differential 
equations ( i.e. Painleve equations) satisfied by the fourth and fifth Painleve 
transcendents are the following : 

p * {z) = + l Pl{z) + AzPl{z) + 2( " 2 " a)p " {z) + R&y (2,15) 

P"(z) - + 1 )P>( Z ? - l -pi( z) + (P^)-l) 2 ( aPi(z) + b 

^ [Z) - [ 2P 5 (z) + P 5 (z)-1 } 5{Z) z n[Z)+ z 2 [ P 5 (z) ] 

(2.16) 

cP 5 (z) dP 5 (z)(P 5 (z) + l) 
z + P 5 (z)-1 ■ 

The Eq.(2.15) and (2.16) are two of the six Painleve equations that serve to 
define the six Painleve transcendents [38]. The fourth and the fifth Painleve 
transcendents allow solutions in terms of rational functions or special func- 
tions [39]. These systems given the Eq.(2.13) and (2.14) are related with 
higher order supersymmetric quantum mechanics (SUSYQM) and we refer 
the reader to the ref.ll, 16, 21 and 22 for more details concerning the ex- 
plicit form of the supercharges, ladder operators, zero modes of the raising 
and lowering operators and the corresponding polynomial Heisenberg alge- 
bras. They possess respectively a polynomial Heisenberg algebra of order 
three and four. The energy of the zero modes are written respectively in 
terms of the parameters of the fourth and fifth Painleve transcendents. The 
corresponding superintegrable systems were discussed in Ref. 14, 16 and 22. 
The Painleve transcendents appear in many context in physics and the re- 
lation with quantum systems with ladder operator and superintegrability is 
very interesting. The study of the classical equivalent of such system is thus 
important and is an unexplored subject. 

In context of classical mechanics, the study of systems with higher order 
ladder operators (i.e. order higher than 2) is a relatively unexplored subject. 
System with third order ladder operators in classical mechanics were explored 
by Eleonsky and Korolev [18], and also in context of superintegrable systems 
[24]. Let us present the system considered in Ref. 24 : 

u 2 

V(x) = —(2b + 5a; 2 + eAxVb + x 2 ). (2.17) 
18 
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A new family of superintegrable systems [21] were also constructed using 
the existence of ladder operators. We will present the general systems possess 
a third order ladder operators and classify systems with fourth order ladder 
operators. As for higher order superintegrability, the classical and quantum 
systems with higher order ladder operators differ. 



3 Classical systems with ladder operators of or- 
der 3 

From Eq.(2.2), a general third order raising and lowering operators has the 
form : 

X I'" ■ <hU)I y ; + h{x)P x ± if (x). (3.1) 

From Eq.(2.5) we obtain the following differential equations : 

-iu + if' 2 (x) = 0, (3.2) 

-uf 2 (x)-f' 1 (x) + 3V'(x) = 0, (3.3) 

- iw/i(x) + if {x) - 2if 2 (x)V'(x) = 0, (3.4) 

-uf o (x) + f 1 (x)V'(x) = 0. (3.5) 

The first two equations give us ( the terms that contain the two integration 
constants can be removed by an appropriate translations of x and V(x) ) : 

f 2 (x) = ux, h{x) = -\j 2 x 2 + 3V(x). (3.6) 

When we put this in the last two equations we obtain : 

- uf' Q + V\x){-u 2 + 3V\x)) + {~u 2 x 2 + 3V{x))V\x) = 0, (3.7) 

/o = -(-wV + 6uV(x) + AuxV\x)). (3.8) 
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We can obtain a second order nonlinear differential equation for V(x) 



l ^-3u 2 V(x)-3u 2 xV'(x) + 3(V'(x)) 2 -^u 2 x 2 V"(x) + 3V(x)V"(x). (3.9) 

This equation coincide with the one obtained in Ref.18 (i.e. the Eq.5.3) 
and the one obtained in context of 2D superintegrability with second and 
third order integrals by Gravel [13] (more precisely the Eq.(19) in the limit 
h — > 0). As shown in Ref.13 and also Ref.18, this equation can be integrated 
by using the fact that it possesses a first integral. 

However instead of giving the solution directly, let us take an other approach 
and impose the constraint that the ladder operators satisfy a polynomial 
Heisenberg algebra. The problem is reduced to solve a system of algebraic 
equations instead of a system of differential equations. Thus we consider the 
second relation of the Eq.(2.6). The polynomial Q{H) can be at most a third 
order polynomial of the Hamiltonian with the following form : 

A X A+ = Q(H) = 8(H - e 1 )(H - e 2 )(H - e 3 ). (3.10) 

Instead of a set of differential equations as the one obtained from Eq.(2.5), 
the Eq.(3.10) generates a set of algebraic relations among fo(x), fi(x), f 2 (x) 
and V(x) involving three parameters ( ei, e 2 and e 3 ) : 

2e 1 + 2e 2 + 2e 3 + 2f 1 (x) + f 2 (x) - QV(x) = 0, (3.11) 

-Ae 1 e 2 -Ae 1 e 3 -Ae 2 e 3 +f 2 (x)+2f (x)f 2 (x)+8e 1 V(x)+8e 2 V(x)+8e ;i V(x)-12V(x) 2 = 0, 

(3.12) 

8eie 2 e 3 + / 2 (x) - 8e l e 2 V(x) - 8e l e z V{x) - 8e 2 e 3 V(x), (3.13) 
+8e 1 V(x) 2 + 8e 2 V(x) 2 + 8e 3 V(x) 2 - 8V(x) 3 = 0. 

We can use the expressions given by Eq.(3.6) for j\ and f 2 and insert them 
in Eq.(3.11), (3.12) and (3.13). We obtain e\ + e 2 + e 3 = and the following 
two equations with 2 parameters c and d ( with d = 4(ef + Ei€ 2 + e 2 ) and 
c = 32u 2 (e 2 e 2 + e l e 2 )) : 
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4^4 



U X 



+ d + 2ooxf (x) - 3uj 2 x 2 V(x) - 3V(x) 2 = 0, (3.14) 



4 

"4^2° + f ^ + 2dV ^ ~ 8V ^ 3 = °" 
The solution is thus (including the Eq.(3.6)): 

fo^ = i-V^) + ^xV(x)-l^-^- x , (3.15) 

-9^ 4 (x)+14w 2 xV 3 (x) + (6d-ywV)y 2 (x) + (^ 6 x 6 -2d)F(x)-(^wV+ ( i 2 ) = 0. 

(3.16) 

This solution satisfy all the equations given by Eq.(3.2)-(3.5) and (3.11)- 
(3.13). This solution agree with results obtained in Ref.18 and 13. The 
Eq.(3.14) coincide with Eq(5.5) of Ref.18 (up to a translation of the coordi- 
nates x and the potential. The Eq.(3.16) is a polynomial equation of degree 
4 and thus can be solved. However, let us do not present explicitly this solu- 
tion and only present a particular case. By imposing two of the roots of the 
polynomial to be equal (i.e. €\ = e 2 ) or ( c 2 = 64 ^" ) we obtain the following 
solutions 

fo = -Qe 2 ux + u 3 x 3 , V(x) = -2e 2 + ^oj 2 x 2 , (3.17) 

and 

J 27 3ux 27 V 

(3.18) 

,2^2 



^>LU T 2 

= 2e 2 + ± -VlSeaX^ + ^x 4 . 

18 9 



We have the following relation e 2 = y j^- We can also use the parameter 

4 7 

d = ^ to recover the potential obtained in Ref.13. 
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4 systems with fourth order ladder operators 



Let us consider the ladder operators of the following form (i.e. as given by 
Eq.(2.2)) : 

= ±iP* + f 3 (x)p 3 ± if 2 (x)p 2 + h{x) V ± if {x). (4.1) 

We obtain from Eq.(2.5) the following five differential equations : 

-u-f 3 {x) = 0, (4.2) 

- iuf 3 (x) + if 2 {x) - UV\x) = 0, (4.3) 

- u>f 2 (x) - f[{x) + 3f 3 (x)V'{x) = 0, (4.4) 

- iuh(x) + if' (x) - 2if 2 V'{x) = 0, (4.5) 

-uf (x) + f 1 (x)V'(x) = 0. (4.6) 

From the Eq.(4.2) and (4.3) we obtain ( again by an appropriate transla- 
tions of x and V(x) we can remove the terms involving the two integration 
constants ): 

h = ~\w 2 x 2 + AV(x), f z {x) = -ux. (4.7) 

We insert the expressions given by Eq.(4.7) in the Eq.(4.4), (4.5) and (4.6). 
From the Eq.(4.6), we can isolate the functions fo(x) and f {x) in terms of 
fi(x), f'i(x), V(x) and V (x). Using these expresions, we obtain from the 
Eq.(4.4) and (4.5) : 

- AuV{x) - f[(x) + ^ux{u 2 x - QV'{x)) = 0, (4.8) 

{u 3 x 2 - 8uV{x) + f x {x))V(x) + fi(x){-u 2 + V"{x)) = 0. (4.9) 

From the Eq.(4.8) and the equation obtained by taking the derivative of the 
Eq.(4.8), we can isolate f[(x) and f[ (x). We can thus insert these expresions 
for fi(x) and f[(x) in the equation obtained by isolating fi(x) in Eq.(4.9) 
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and taking the derivative. We thus obtain for V(x) the following third order 
nonlinear differential equation : 

oo 4 x 2 +2(-4u 2 V(x)-Qu 2 xV' \x)+(15V' ' (x)) 2 +(16V(x)+x(-2oj 2 x+9V' \x)))V \x)) 

(4.10) 

3V'(x)(8V(x) + x{-u 2 x + 2V'{x)))V"'{x) 
+ u 2 - V"{x) ~ ' 

with u 2 - V"(x) ^ 0. The solution V" (x) = oj 2 is the well know harmonic 
oscillator. 

Again, let us take the same approach than for systems with third order lad- 
der operators and impose for the raising and lowering operators to generate 
a polynomial Heisenberg algebra (i.e. impose the Eq.(2.6) to be satisfied). 
The product of these opperators can be at most a quartic polynomial of the 
following form : 

A X A + X = 16(H - ei )(H - e 2 )(H - e 3 )(H - e 4 ). (4.11) 

Using the explicit form of the ladder operators given by Eq.(4.1), we obtain 
the following set of equations : 

2( ei + e 2 + e 3 + e 4 ) + 2/ 2 (x) + / 3 2 (x) - 8V(x) = 0, (4.12) 

- 4(eie 2 + eie 3 + e 2 e 3 + eie 4 + + e 3 e 4 ) (4.13) 
+2/ (x) + f 2 (x) + 2/ 1 (x)/ 3 (x) + 12( Cl + e 2 + e 3 + e 4 )V(x) - 24V(x) 2 = 0, 
8(eie 2 e 3 + ei e 2 e 4 + ei e 3 e 4 + e 2 e 3 e 4 ) + f 2 (x) + 2f (x)f 2 (x) - 16(eie 2 + e x e 2 (4. 14) 
+e 1 e 3 + e 2 e 3 + e 1 e 4 + e 2 e 4 + e 3 e 4 )U(x) + 24(e 1 + e 2 + e 3 + e 4 )U(x) 2 -32y(x) 3 = 0, 
- 16eie 2 e 3 e 4 + f (x) 2 + 16(eie 2 e 3 + eie 2 e 4 + eie 3 e 4 + e 2 e 3 e 4 )V(x) (4.15) 
-16(e 1 e 2 +e 1 e 3 +e 2 e 3 +e 1 e 4 +e 2 e 4 +e 3 e 4 )y(x) 2 +16(e 1 +e 2 +e 3 +e 4 )\/(x) 3 -16y(a;) 4 = 0. 

We will insert the expression given in Eq.(4.7) for / 3 (x) and / 2 (x) in 
Eq. (4. 12)- (4. 15). We obtain the following relation e\ + e 2 + e 3 + e 4 = 
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and using the following parameters : 

d=e\ + ei e 2 + e 2 + ei e 3 + e 2 e 3 + e 2 , (4.16) 

c = e 2 e 2 - e t 4 - 4 e 3 ~ 2eie 2 e 3 - e 2 e 3 - €163 - e 2 e 2 , (4.17) 

e = eie 2 e 3 + eie 2 e 3 + eie 2 e^, (4.18) 

we obtain three equations for fo(x), fi(x) and V(x) : 
4 4 

- d + + 2/ (x) - 2wx/i(o;) - 4u 2 x 2 V{x) - 8V 2 {x) = 0, (4.19) 

8c - u 2 x 2 f (x) + f 2 (x) - AdV{x) + 8f (x)V(x) - 32V(x) 3 = 0, (4.20) 

- 16e + / (x) 2 + IQcV(x) - 4ciy(x) 2 - lQV 4 (x) = 0. (4.21) 

The solution is thus (including Eq.(4.7)) : 



f (x) = y/-16e + 16cV(x) - 16dV(x) 2 + 16V(x) A , (4.22) 

f^x) = — (I6d + uV + 8/ (x) - 16u 2 x 2 V(x) - 32V 2 (x)), (4.23) 

(4.24) 

V(x) = -32 W 2 xV 5 (x) + (128d+17wV)y 4 (x) + (-96 C -^^-^-32^V-^^)\/ 3 (x) 

/ - ,9 256c 2 256cd 9 9 ,44 llu^x®,.,. s 

+ -40d 2 + 352e + — - + + 32cu 2 x 2 + 3du 4 x 4 + oo )V 2 {x) 

urx 4 urar 32 

, 128cd 2 512ce 256c 2 64d 3 256de co w x 10 ^ r/ , 
+(-16cd ri j-j- + -3-2 + —2-3- 



ura; 4 w 4 ar o> 2 a; 2 cj 2 ar urar 64 
16d 4 128d 2 e 256e 2 64cci 2 256ce 



+64c 2 + 4cf - 112de + + - " + -^j- - : ^ rv + - 8cdu 2 x 2 

o> 4 x 4 u; 4 ar w 4 x 4 urar urar 

+ ^2 4 4 17 4 4 _ 1 6 X 6 + 1^8 ^! 

2 2 4 64 4096 

These relations satisfy well the Eq.(4.2)-(4.6) and Eq.(4.12)-4.16). Let us 
again do not present explicitly the solution for the potential that can be 
obtained because this is a quintic equation. We present a particular solution 
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by considering a triple roots (i.e. €\ = e 2 and e 3 = e 2 or c = |y |d2 and 

C 12 ^ 

We obtain the following solution : 

(-8e 2 + a;V) (-8e 2 + u,V)3((8e 2 + u,V) 
4^ ' /o(x) - 16^ ' (425) 

V(x) = -e 2 + ^wV + 

and 

the following ID system : 

3 7 1 

/ (x) = -3e^+e 2 (--w; 2 x 2 ±— ^Q4e 2 u 2 x 2 + u 4 x 4 ) + — (-15cu 4 x 4 +17u 2 x 2 ^6Ae 2 u 2 x 2 + u A x A ) 
2 16 512 

(4.26) 

f(\ r 3w3x3 + e 2 V646 2 q;^ + a; 4 x 4 5 /- 2 4 4 

ji(x) = -6e 2 ux — — ± ± — ux\/<oAe 2 u 2 x I + u 4 x 4 , 

32 ojx 32 

(4.27) 

V(x) = — (96e 2 + 5w 2 x 2 =f 3 ^64e 2 oj 2 x 2 + u 4 x 4 ). (4.28) 
64 



5 New superintegrable systems in classical me- 
chanics 

Let us now consider a classical 2D Hamiltonian separable in Cartesian coor- 
dinates 

H(x 1 ,x 2 ,P 1 ,P 2 ) = H 1 (x 1 ,P 1 ) + H 2 (x 2 ,P 2 ), (5.1) 

for which polynomial ladder operators (A x . and A£ ) of order rij exist and 
satisfy polynomial Heisenberg algebras. From the Eq.(2.2) the operators 
fi = (A+) mi (A") m2 and f 2 = (A") mi (A+) m2 ( of order n x m x + n 2 m 2 ) 
Poisson commute with the Hamiltonian H given by Eq.(2.1) if m,\U xl = 
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^2^x2 = w with mi ,m2 £ The following products are also polynomial 
(i.e. polynomial in the momenta) integrals of the Hamiltonian H : 

h = «) mi (^ a r-(^~ r«r, / 2 = «) mi (^) m9 +(^) mi (^i) ma • 

(5.2) 

By construction the Hamiltonian has the following second order integral re- 
lated to separation of variables in Cartesian coordinates : 

K = H 1 - H 2 . (5.3) 

The Hamiltonian H given by Eq.(5.1) is thus maximally superintegrable (i.e. 
there are 3 functionally independent integrals of motion including the Hamil- 
tonian himself). We construct from integrals given by Eq.(5.2) and (5.3) the 
following polynomial Poisson algebra : 

{K, h} p = 2ul 2 , {K, I 2 } p = 2uh, {h, I 2 } p = 2Qi(i(# + K)) m ^ 1 

(5.4) 

5.1 Two new families of superintegrable systems 

We can consider a new following 2D superintegrable system H as given by 
Eq.(5.1) with H\ and H 2 given by : 

p2 p2 

Hi = ^ + V 1 (x 1 ; wi), H 2 = + V 2 (x 2 ; cu 2 ), (5.5) 

with Vi(xi,Ui) and V 2 (x 2 ; u 2 ) that both satisfy Eq.(3.16). We impose the 
constraint miU) Xl — m 2 u X2 = with mi,m 2 G Z + . From the construction 
at the beginning of the Section 5, this systems is thus maximally superinte- 
grable. 

Using the same form (i.e. the one given by Eq.(5.5)), we impose that V\ 
and V 2 satisfy now the Eq.(4.24). We can thus generate a second family of 2D 
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Figure 1: 



Figure 2: 



superintegrable systems. These results show how new systems with ladder 
operators can generate new superintegrable Hamiltonian. These results can 
be extended in N dimensions. 

Let us present trajectories for a particular case (i.e. satisfying Eq.(4.28)) 
respectively for e = 1 and e = — 1. All bounded trajectories are closed and the 
motion is periodic. We choose for Fig.l and Fig.2 the following parameters 
wi = l,u 2 = 2, e 2xi = 4, e 2x2 = 16, x(0) = y(0) = 1, x'(0) = 1, y'(0) = -3 
and t = [0,20]. 

6 Conclusion 

We presented the most general ID classical system with a third order ladder 
operators and also the most general system with fourth order ladder opera- 
tors. We obtain respectively potentials that satisfy a fourth and a fifth order 
algebraic equations. They are the classical equivalent of the quantum sys- 
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terns written in terms of the fourth and fifth Painleve transcendents obtained 
and studied in earlier papers [11,16,19,21,22]. This point out how systems 
with higher order ladder operators differ in classical and quantum mechanics. 

Using a method introduced in earlier articles [14,24] and extended by Miller, 
Kalnins and Kress [36] to systems with separation of variables in polar co- 
ordinates, we obtained two new families of superintegrable systems that are 
deformations of the anisotropic harmonic oscillator and Evans- Verrier sys- 
tems [40]. This paper point out also how ladder operators can be used to 
generate new superintegrable systems. Let us mention other recent papers 
that discuss ladder operators in context of superintegrability [41,42]. 
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